International Journal of Engineering and Applied Sciences (IJEAS) 

ISSN: 2394-3661, Volume-2, Issue-12, December 2015 


Multivariable Control of a Binary Distillation 

Column 

Jose Leonardo Benavides Maldonado, Jorge Michael Valarezo, Darwin Tapia Peralta 


Abstract — this work is present students of University National 
of Loja in the southern Ecuador and associated people with the oil 
industry. Two strategies to control some parameters of a binary 
distillation tower, such as LQR (Linear Quadratic Regulator) and 
another in intelligent control such as fuzzy logic controllers 

With purpose of comparison with the help of MATLAB® 
software and Simulink® the benefits offered by the use of either 
strategy in the oil industry. 

Index Terms — Chemical Industry, Distillation Columns, 
Fuzzy Control, LQR Control. 

I. Introduction 

This The need for professionals who can work in the oil 
industry, is growing in Ecuador so it is important that students 
of the National University of Loja will relate more to this 
process in their academic training so one of the objectives of 
this paper is intended to raise awareness of this process. 

An application to understand this process is a simple 
classical binary distillation tower where unlike those used in 
the oil industry is responsible for separating a liquid mixture 
of two substances with different degrees of volatility 
document is a template. An electronic copy can downloaded 
from the conference website. For questions on paper 
guidelines, please contact the conference publications 
committee as indicated on the conference website. 
Information about final paper submission is available from 
the conference website. 

This achieved because a distillation tower is a system for 
continuous reheating reheated at the bottom [1]. This shown 
more detail in Fig. 1. Where it shows a distillation column in 
which a flow F with Zi composition is fed, a distillate D in the 
top you will have an Xd composition and B flow with Xb 
composition also can be seen a reflux L. 

This will always exist for a more adequate purification, this 
reflux is the part that is condensed and the other is the part that 
comes out as a finished product this part of the tower. 

Known as rectification zone here increases the more 
volatile component, while at the bottom the heaviest component 
increases and thus reduces the lighter component [2]. 

Can also see in this same figure on the right a reflux drum 
where the condensate temporarily stored, to distill and a pump 
to evacuate the above [3]. 
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Fig. 1 In the left figure the stages shown operating of a binary distillation 
tower and in the right of this figure a distillation tower with a drum to control 
the level of reflux. Source photograph: left 

http://www.voutube.com/watch7vHeCZ4UKfLKc , right 

http://www.youtube.com/watch?v=sYIkvDOVTzA&list=PL995A30D4624 

4C35B 

II. Materials and Methods 

For nonlinear systems there are a large variety of structures 
and techniques such as, neural networks, support vector 
machines least squares (LVS-SVMs), fuzzy logic [4]. Two of 
these strategies have elected to meet certain parameters of the 
distillation tower. 

2.7 Definition of state space 

The idea comes from paper the state-space method for 
describing in state variable differential equations. Here the 
dynamical systems described by a set of differential equations 
of first order in known state variables and the solution can 
considered as a solution in the state trajectory. This method is 
ideal to perform calculations with the help of a computer. The 
design of state-space control is a technique by which the 
control engineer designs the compensation directly using the 
description of the system state variables. 

Also called modern control design, while methods based 
on transfer function, as the locus of the roots and the 
frequency response called classical control design, but their 
release date is a mistake to refer as modern, and it is 
preferable to refer to the two as an approach to the design and 
methods of state space method [5]. 

The stable models are described by two equations the state 
equation which is shown by one, two. 

X=Ax+Bu (1) 

y = Cx + Du (2) 

Where Y (t) represents the output of the system, x (t) is the 
system state, u (t) = the system input and X is the future state 
of the system depending on the current state and of the output 
of the input of the system [6]. 
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2.2 Linear Quadratic Regulator (LQR) 

Assuming all the system state can measured, we know 
possible to design a state feedback control that minimizes the 
value of a quadratic cost functional infinite horizon defined 
by three. 

00 t 1 T ( 3 ) 

J := J ( X Qx+ e u u)dt ' ' 

0 

Where Q is a positive matrix and is a positive constant, the 
control law shown in equation four. 


u _opt = —r 1 B t Px (4) 

Where P is positive definitive matrix as a solution to the 
algebraic Riccati equation (ARE) with equation five. 

a t p + pa-^pbb t p-q = o (5) 



Fig. 2 Binary distillation column with pressure variation 

The equation of binary distillation column with n plates for 
a general binary system of components given in the above 
reference six, seven, eight. This following linearized model 
obtained from the above reference, for a column containing 
eight plates 

X = Ax+Bu+Ew ( 6 ) 


That can be proved to exist as soon as A, B is a stabilizable 
pair, the firsts and seconds terms correspond to the energy of 
the controlled output and control signal, respectively, and the 
LQR strategy is to minimize , however, the decrease of a 
requires of both. The paper is to establish a dot between these 
behaviors in the conflict that given by the value of e [7]. 

2.2 Fuzzy logic-based controllers 

Fuzzy logic aims to give a greater range of options when 
finding a solution to a particular problem but the solution 
given to a problem does not always allow a greater horizon of 
probabilities, 

Zadeh proposed this problem the resolver by applying the 
logic of Lukasiewicz to find a solution designating each 
object in a set and I think a complete algebra for fuzzy sets. 

This theory proposed membership functions (p), or false 
and true values over the range [0.0, 1.0], the fuzzy degrees are 
not the same as the odds, these measure indicate something 
will happen or not, while fuzzy measure levels are a condition, 
i.e. has to be given a degree of belonging to the elements. The 
second step is to propose linguistic rules (inference), the 
degree of membership of each of the variables is evaluated on 
a set of inference rules. These were determined with the help 
of an expert. The inference rule set determines a consequence, 
that is, assigned a degree of membership of each fuzzy set, 
which characterizes the output. 

After obtaining the consequences, the third stage is a 
process to determine the optimal values of output, known as 
defuzzification, 

Which involves passing the grade of membership, from the 
result of the inference rule to a crisp or actual value to do this, 
previous membership functions of each of the outputs tuned in 
order to obtain a measurable value [8]. 

2.3 Description of the process 


y = Cx 



ym = x 



Where u, is the input: 

u=ul, u2, u3, (ul= Reboiler steam temperatura; u2= 
Condenser coolant temperatura; u3= Controlled reflux). 

Y is the output: 

y=yl, y2, y3, (yl Composition of more volatile component in 
the bottom product; y2=Composition of more volatile 
component in top product; y3=Pressure). 

XI Composition of more volatile component in condenser 
XIOComposition of more volatile component in reboiler 
XI l=Pressure 

X2Composition of more volatile component in the plate #1 


X9=Composition of more volatile component in the plate #8 
W1 Change of input feed concentration 
It desired to design a controller to regulate three outputs yl, 
y2, y3 against the unmeasurable disturbance W1 and other 
unmeasurable disturbances, with as fact a settling time as 
possible, subject to the constraint. 


\u x \ <2.5 


Ur 


U Z 


<2.5 

<0.30 


t> 0 


For W 1=1 in this problem, the choice of what measurable 
outputs to use in the controller, is considered to be part of the 
problem statement; in general, a controller which uses the 
fewest number of output measurements is desirable [9] . 


2.4 Process variables 

The process under study is a distillation tower whose 
variables are; Temperature, Pressure, Level and Flow, 
variables involved in the analysis are input (ul, u2, u3). 


The problem describes a realistic problem of a binary 
distillation column, and has the feature that pressure variation 
is included in the model’s description: the system is 
multivariable, with three inputs and three outputs, and 
includes one disturbance input; this presented in Fig. 2. 


Which are replaced by the following abbreviations (Tf, Tr, 
Re) and output (yl, y2, y3) which are replaced by the 
following abbreviations for greater understanding (Cf, Ct 0 , 
Pt), which can be seen in Table 1. 
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Table 1: Process Variables 


MODELO DE LA PLANTA 

Input variables 

Output variables 

Reboiler steam temperature ( Tf 
) [ °C] 

Composition of more volatile 
component in the bottom product 
( Cf) [ g/ml xlOO =%, or mole 
fraction is dimensionless] 

Condenser coolant temperature ( 
Tr ) [°C] 

Composition of more volatile 
component in top product 
( Cto ) [g/ml xlOO =%, or mole 
fraction is dimensionless] 

Controlled reflux (Re) [ g/min] 

Pressure in the top (Pt) [ which is 
the same or Pascal Pa] 

Disturbances 

Change of input feed 

concentration 
( W7 ) [g/min] 


Then in in Fig. 3. The system variables represented in a 
black box model where it is clearly a model is MIMO or 
multivariable system with three inputs, three outputs and one 
disturbance that is food, the mixtures is alcohol and water. 



Fig. 3 Black box model of Binary distillation column, Source photograph: 
left [9], right [4] 

When to wear this in state space system and simulate in 
Simulink, to understand the behavior of the outputs (Cf) 
(CTO), (Pt) to the inputs (Tf-2.5), (Tr-1.5), (Re-0.30) type 
step, is that the system is unstable as shown in Fig. 4. 



Fig. 4. Characteristics of the output of the Binary distillation column 


To correct this unstable behavior of the system shown in 
the figure above, should be developed a small program in the 
MATLAB® editor that allows to present this state-space 
representation to an equation transfer function obtained will 
be carried Simulink® to be controlled through the design of 
three fuzzy controllers and then is one of the objectives of this 
work. 

III. Results and Discussion 
3.7 Control based on LQR 

Is carried out a program in MATLAB® editor, using the 
toolbox with predictive models that have, specifically the ss 
command to state space, Here 4 matrices of different order are 
defined (Matrix A, B, C, D, and a sampling period). 

These regulators LQR they have a penalty of exit and entry 
by the R matrix, and outputs the matrix u, and then a sequence 
of operations that must be performed comes, 

Here it is observed that the polynomial of Y has a penalty 
Q, and the polynomial of u has an R penalty and then this 
multiplied by Q and that is the function to be minimized when 
the minimum of this function is optimal [10], [11]. 

The Fig. 5. Gives the result of LQR strategy applied to the 
process, indicating that all three outputs converge to zero as 
expected except that y3 is minimal but the difference to zero 
is6e~ 4 [12] [13]. 



Fig. 5 Response to the three outputs of the system and as you can see is stable 
because all the answers emerge to zero 

The response to the statements with a quadratic control 
strategy is almost zero; this can see in Fig. 6. 



Fig. 6 Response of the 11 intermediate states when are use a LQR strategy of 
control 
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3.2 Control based on Fuzzy Logic 


Table 2: Fuzzy rules applied to variables input and output of 

the process. 


For controller design Fuzzy must take into account some 
considerations in the distillation vapor in the process, such 
that the steam flow can considered constant throughout the 
column, 

It is considered the evaporation process on each plate, 
overflow equimolar, which means that each molecule 
condenses delivery equal amount of energy which it absorbs a 
molecule liquid evaporates, thus the seam flow leaving the 
re-boiled, it is the same to the condenser, though with 
different concentrations through the variable control system 
shown in Fig. 7. 



D 

X. 


./a 

titi 


Fig. 7 Control scheme in where shown the sensors and actuators that must 
handle with the Fuzzy controller for output values desired. Photo 
Source: http://www.youtube.com/watch?v=sYIkyDQVTzA&list=PL995A3 
0D46244C35B 

3.2.1 Fuzzy Logic in MATLAB® and Simulink® 

For the development of the rules of inference were used the 
criteria found in the previous tutorial [14], as well as data 
obtained from [4]. 

As mentioned above, the modeling has raised out using 
fuzzy logic techniques, due to their ability to model the human 
procedural knowledge and has been prepared by supporting 
the tool MATLAB®. Work has done using a type of Mandani 
driver because it provides great advantages [15]. 

Then with the help of MATLAB® and your toolbox fuzzy 
you choose the functions membership and fuzzy rules are 
these can be seen in Table 2. 


Variable 

Label Language 

Reflux(0-2.5) 

ConCooTem( 

0-2.5) 

ReboSteam 

Tem 

(0-0.30) 

Input 

Feed(O-l) 

Comdes(0 

-0.08) 

LITTLE 

HIGH 

LITTLE 

HIGH 

LITTLE 

HIGH 

LITTLE 

HIGH 

LITTLE 

HIGH 

Comsed(0 

-0.02 

HIGH 

LITTLE 

LITTLE 

HIGH 

LITTLE 

LITTLE 

HIGH 

HIGH 

HIGH 

HIGH 

Pre(O-l) 

HIGH 

HIGH 

LITTLE 

HIGH 

LITTLE 

LITTLE 

LITTLE 

HIGH 

LITTLE 

HIGH 


In table 3 shows the results of the fuzzy rules of the process 
that must comply with the principle that the distillate is a very 
high concentration approximately 1 being the pure the most 
volatile element pure and the pellet has a concentration near 
shown so zero less volatile component. 

Table 3: Simulated test for results 


Variabl 

e 

Reflujo(0- 

2.5) 

Temsalco 

n(0-2.5) 

Temreher 

(0-0.30) 

Flujodeen 

(0-1) 

Comdes(0 

-0.8) 

Corns 

ed(0- 

0.02) 

Pre(0- 

1) 

Comde 

s(0-0.8) 

0.43 

2.02 

0.0121 

0.359 

0.887 



Comse 

d(0-0.0 

2 

0.161 

0.193 

0.242 

0.0154 


0.015 

1 


Pre(0-1 

) 

1.9 

0.1 

0.27 

0.1 



0.792 

Pre(0-1 

) 

0.305 

0.909 

0.0158 

0.959 



0.163 


In Fig. 8 are shown the subsystem that was conducted to 
obtain the results of the figures 10, 11 y 12. 


irtiMAImmilWil 


nH 


IF 



□ 




Fig. 8. Process represented as a subsystem, applying MATLAB 


® 


In the Fig. 9, the control development with the help of 


Simulik® schown to control nine functions of transfer 
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Fig. 9. Control Subsystem using Simulink® 


In the Fig. 10 shown the result the less volatile 
concentration the bottom sediment as close to zero, for this 
case would be the water 
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Fig. 10. Here is showing the concentration less volatile (water). 

In the Fig. 11 shown the control result of the volatile 
composition in the top, the distillate is the concentration more 
volatile with a value close to one 



Fig. 11. Here is showing the concentration more volatile (alcohol). 


In the Fig. 12 shown the result of pressure at the top of the 
column of distillation. 



Fig. 12. Pressure in the top of distillation column binary 
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IV .Conclusion 

1. Are 9 transfer functions are obtained using state space 

the method because this process has three inputs and 
three outputs 

2. The degree of the denominator is equal to the three 

transfer functions and represents the input state space 
or, as has 11 intermediate states X (is of the 11 order) 

3. Because states are intermediate variables, combination y 

= Cx + Du is giving control outputs in this case are 
transfer functions (the nine) are shown in the 
subsystem of figure 8 

4. The use of LQR strategy an allowed to converge to zero 

all the variables to control because it minimized to 
zero the square error (reference least output all 
squared). 

5. The use of Fuzzy controller allows students of UNL had 

better understand the dynamics of distillation column. 

6. Considering all the working conditions of the distillation 

column can assumed that the enthalpy between the 
dishes are equal. 

7. Both methods are very effective in the control of the 

parameters of binary distillation tower but the 
potential of fuzzy method (fuzzy) allows better 
results in nonlinear processes (is presented as 
example the original model of the column.). 

8. The errors between the two Controllers the pressure is 

low in the LQR (0, 0005), while in ours fuzzy 
controllers is lower. 
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